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CN ■ Abstract. We summarize new pQCD results on pion production in proton- 

proton (pp), proton- nucleus (pA) and nucleus-nucleus (AA) collisions. Our 
calculation introduces intrinsic parton transverse momentum (hr) and is per- 
formed effectively at next-to- leading order (NLO), applying a K factor ex- 
' q ■ tracted for jet events. Two different factorization scales, Q = PT.jet/2 and 

^ | Pt jet are used. Experimental data in pA collisions imply a preference for the 

Ch ■ latter choice at NLO level. We display our results at CERN SPS for AA colli- 
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1. Introduction 

Recent RHIC and CERN experiments require solid theoretical baseline calculations, 
able to reproduce proton-proton, proton-nucleus and the measured nucleus- nucleus 
data accurately enough to find possible deviations due to new collective phenomena 
or to a new state of the matter. In this paper pion production is calculated in a 
pQCD-improved parton model [ |l[. Recently we have performed a leading order 
(LO) calculation and displayed the results of a systematic analysis of pp (pp), pA 
and AA data on pion production [ || . The key points were the introduction of the 
intrinsic transverse momentum distribution of the partons and the folding of nuclear 
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multiscattering effects into the pQCD calculations. We have also calculated pion 
production in next-to-leading order (NLO) [ [|. However, the NLO cross sections 
were estimated by applying an oversimplified constant K factor (K = 2) to the 
LO cross sections. Since the determination of nuclear effects (e.g. the strength of 
nuclear multiscattering) and other collective effects require great precision (better 
than 10 — 20 %), we need to repeat the earlier calculation using appropriate NLO 
cross sections. In Section 2 we summarize our results on energy and momentum 
dependent K factors. 

In standard calculations for pp collisions the factorization scale is typically cho- 
sen in the range Pt/3 < Q < 2pT- We find that at the NLO level the reproduction 
of nuclear effects in pA collisions (specifically the Cronin effect [ |2| ||) strongly 
depends on the choice of the factorization scale. We discuss this issue in Section 3. 




Pw (GeV) 

Fig. 1. The Kj et (s,pTj e t) factor as a function of yfs and PT,jet a t scale Q = 
PT,jet/2- Solid lines guide the eye through the calculated points and the fitted 
parabolae of eq. (g) are indicated by dashed lines [ |) . 



2. K factor for NLO jet-production in pp collision 

In a previous paper [ |6| we have extracted a K factor for jet production in pp 
collisions. The scale in the parton distribution functions (PDFs) was chosen to be 
Q = pxjet/2 for jets with jet-cone angle R = 1 and separation R sep — 2R. An 
energy and transverse momentum dependent correction factor, Kj et (s,pTj e t), was 
then obtained as 

da NLO dcr Born 

d £ = K jet (s,PT,jet) -g , (1) 

where the cross sections are understood at the jet level (see eq. (S)). 
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Full NLO calculations were performed by the Ellis - Kunszt - Soper (EKS) [ [?], || 
and the Aversa-Chiappetta- Greco -Guillet (ACGG) [ |j| groups, and are based 
on the matrix elements published in Ref . [ |l(| . We used a public Fortran code from 
the EKS group [ with MRST (central gluon) PDFs [ pi to calculate the Born 
and the NLO contributions. The extracted K factor of eq. (Ill) is shown in Fig. 1 in 
the energy range 20 GeV < y/s < 200 GeV. 

For practical purposes, the K factors were parameterized (within 2-4 % preci- 
sion) as 

K jet (s,p T , jet ) = 1.6 + H| _ _^__ PTjet + __j__ p | , ei , (2 ) 

where Prjet and ^fs are in GeV and the constants are understood with their ap- 
propriate units. 

In this paper we plan to use a different scale in the PDFs (see Section 3), namely 
Q = pT.jet- Thus we repeat our calculation of the jet-level K factor at this new 
scale, using R = 1. The results are displayed in Fig. 2. 




1 00 GeV 

at:""" 140 GeV 



200 GeV 



n Pw (GeV) 

Fig. 2. Same as Fig. |l|, for scale Q — pr.jet- 

These new results may be again parameterized quadratically with good precision: 

20 45 T 

K jet (s,Pr, jet ) = 1.79+ - {V _J 7 )2 p Tjet + (V - s : Q )2 phet • (3) 

Figures [j] and || demonstrate that increasing the scale by a factor of 2, the 
jet-level K factor is increased by ~ 10 — 15 %, maintaining the characteristic 
dependence on s and pt- At RHIC energies (-^/s = 200 AGeV) the K factor is 
almost constant with values ~ 1.7 (at Q = prjet/2) and ~ 1.9 (at Q = prjet)- 
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3. Pion production from CERN to Tevatron energies 

Here, we summarize the method used to calculate pion production in a pQCD- 
improved parton model, following our previous work in LO, see Ref. [§. In Subsec- 



tion 3.1 we introduce the intrinsic transverse momentum distribution for partons 
and fit the width of the applied Gaussian distribution to available pp (or pp) data 
in the cm. energy range 20 GeV < y/s < 1.8 TeV. A large collection of pion and 
charged hadron production data can be utilized in this energy range to extract 
information about the intrinsic transverse moment um distribution, e.g the energy 



dependence of the Gaussian width. In Subsection 3.2 the Cronin effect [ f|, |) is 



discussed in pA collisions. Subsection |3.3| deals with hard pion production in AA 
collisions at CERN SPS energies. 

3.1. Parton model calculations with intrinsic transverse momentum 

In pp collision we describe the invariant cross section for pion production in a 
pQCD-improved parton model on the basis of the factorization theorem [Qj: 

E TT—^r- = dx a dx b dz c f a / p (x a ,Q 2 ) f b /p(x b ,Q 2 ) x 

d P tZTJ 

D„ /c {z c ,Q' 2 ) 



abed L 



ab^cd 

K Je t(s,PT,jet) -TJ 

at 



s S(s + t + u) . (4) 



Here f a /p(x a ,Q 2 ) and fb/ P (xb,Q 2 ) are the LO/NLO parton distribution functions 
for the colliding partons a and b within the interacting protons as functions of 
momentum fraction x, at scale Q. da/dt is the hard scattering cross section of 
the partonic subprocess ab — > cd in LO (Born term) and the Kj et factor is applied 
to include NLO contributions. In eq. (Q) we use the convention that the parton- 
level Mandelstam variables are indicated by a 'hat' (like t above). The LO/NLO 
fragmentation function (FF), D^/ c (z c ,Q' 2 ) gives the probability for parton c to 
fragment into tt with momentum fraction z c and fragmentation scale Q' . For the 
results presented here we fix Q' = pr/2 and use the KKP parametrization [ |l3[]. 

For LO calculations, Kj et (s,pT,jet) = 1- As discussed in Section 2, in NLO 
calculations the Born partonic cross section is multiplied by Kj et (s,pTjet) inside 
the integral of eq. (|). For LO calculation the GRV [ [l4|], for NLO the MRST [ Q 
PDF sets were applied at scales Q = PT,jet/% an d Q — PTjeti respectively, where 

PTJet = Pt/z c - 

In a phenomenological approach, eq. (|]) can be generalized to incorporate in- 
trinsic transverse momentum by using a product assumption and extending each 
integral over the parton distribution functions to fc^-space [ |l5|, [l^] , 

da; f a /p(x, Q 2 ) -> dx d 2 k T g(k T ) f a / p (x, Q 2 ) , (5) 

where g(kr) is the intrinsic transverse momentum distribution of the relevant parton 
in the proton. We follow this approach in the present work, choosing g(kr) to be a 
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s ,/2 (GeV) 

Fig. 3. Best fit values of (k T ) in pp — > ttX or pp — > h^X reactions at LO ( dashed 
line) and NLO (solid line) levels with factorization scale Q = pr,jet/2, both with 
errors (dotted lines). Due to the overlap at the same energy, some points with error 
bars have been slightly shifted for better visibility. Thin lines represent the error 
estimate of the curves, while the thick bars are for RHIC energies in NLO. 



Gaussian: 

9 <« ■ say*** ' (6) 

with (k T ) being the 2-dimensional width of the kx distribution. Increasing the 
value of parameter (fey), the particle production also increases in the transverse 
momentum window 2 GeV < px < 6 GeV. 

In Figs. H and 0, we present the best fit values of (k T ), calculated in pion and 
charged hadron (h ) production in pp and pp collisions at scales Q = Vrjet/2 and 
Q = PTJet , respectively, using several independent experimental data [ [|, [Tt], [l^, 
po| , |22] , p3| , pl| |2^| . As expected, the obtained value of (fcf.) at NLO level is less 
than at LO level if Q — pT,jet/2 is used in both calculations (see Fig. 3). However, 
this difference in (k T ) almost disappears changing the NLO scale to Q = PTjet, but 
keeping the LO scale at Q — pTjet/2 (see Fig. 4). 

As an example to illustrate the degree of accuracy of the description in the pp 
sector, Fig. || compares calculated tt + and tt~ spectra and tt~ /ir + ratios to the data [ 
U at the cm. energy y/s — 27 A GeV, respectively, obtained with the values of (k T ) 
(the different scales are indicated in the top panel). We find that the data/theory 
and 7T~ /ir + ratios are well reproduced for 2 GeV < pt < 6 GeV. Based on this and 
similar examples, we believe that hard pion production in pp collisions is reasonably 
under control at the present level of calculation. 
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3.2. Hard pions from pA collisions - Cronin effect 

Interest in the nuclear dependence of hard particle production was motivated by 
the discovery of the Cronin effect [ ||, |j| : it was found experimentally that in the 
transverse momentum window 2 GeV < p? < 6 GeV more particles are produced 
than it was expected from a simple scaling of the pp data. As discussed in the 
introduction, reproducing the Cronin effect could be the key point - at least in the 
sense of scale fixing - in our calculations. 

In our model there is an extra contribution to the basic Gaussian width of the 
intrinsic parton transverse momentum distribution due to the nuclear environment. 
This extra width can be related to the number of nucleon-nucleon (NN) collisions 
in the medium. To characterize the (k T ) enhancement, we write the modified width 
of the parton transverse momentum distribution of the incoming proton as 

(k T ) pA = {k 2 T ) PP + C- h pA (b) . (7) 

Here (k^} pp is the width o f th e transverse momentum distribution of partons in pp 
collisions from Subsection [O] (also denoted simply by (&§■)), h p A{b) describes the 
number of effective NN collision at impact parameter b, each of which imparts an 
average increase in the transverse momentum width denoted by C . In pA reactions, 
one parton from the projectile proton and another one from the target nucleus 
participate in the hard collision. We apply (k^) p A for the projectile parton, which 
incorporates the additional NN collisions indicated in eq. (f?). For the target parton 
the original {k^) pp is used from Fig. || or Fig. ||. 

The effectivity function h p A (b) can be written in terms of the number of col- 
lisions suffered by the incoming proton in the target nucleus, uaQ)) — <^NNtA{b) 1 
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Fig. 5. LO and NLO calculations at different Scales' (see text). Top panel: 
invariant cross section of 7r + and tt~ production from pp collisions (Negative pion 
data are divided by 10 for better visibility). Central panel: data/theory ratios, 
D/T. Lower panel: ir~ /ir + ratios as functions of transverse momentum at cm. 
energies y/s =27 .4 GeV. Data are from [ p|. 



where unn is the inelastic nucleon-nucleon cross section: 

va{V) - 1 va{V) < V-n 



h P A{b) = \ „ _ 1 otherwise • (8) 



00 corresponds to the case where all possible semihard colli- 
sions contribute to the broadening. Requiring independence of target we found that 
for realistic nuclei the maximum number of semihard collisions is 2 < v max < 5. 

According to the Glauber picture, the hard pion production cross section from 
pA reactions can be written as an integral over impact parameter b: 

e*£ = u% tA{ b) e, ^imuiMM , (9) 

dp J dp 

where the proton-proton cross section on the right hand side represents the cross 
section from eq. (Q) with the transverse momentum extension as given by eq.-s (||) 
and (^). Here = / dzp(b,z) is the nuclear thickness function (in terms of 

the density distribution p) normalized as J d 2 btA(b) = A. Furthermore, the PDFs 
are modified in the nuclear environment ( "shadowing" ) [ [2?], p8| . This effect and 
isospin asymmetry are taken into account on average using a scale independent 
parameterization of the shadowing function S a / A {%) adopted from Ref. [ |2?J : 



fa/A(x, Q 2 ) = S a / A (x) 



A 



fa/ P (x,Q 2 )+ (l-J^j fa/n(x,Q 2 ) 



(10) 
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Fig. 6. Normalized W/Be charged pion cross section ratios at different energies: 
yfs = 19.4; 23.8; 27.4; and 38.8 GeV ||, || ||. The calculations are at NLO 



level with scale Q — PT,jet/2 (dashed lines) and with Q 
deviation from unity represents the Cronin enhancement. 



PTjet (solid line). The 



where f a / n ( x , Q 2 ) is the PDF for the neutron and Z is the number of protons. 

The Cronin enhancement is presented in Fig. ^| by the normalized (with mass 
number ^4) W/Be charged pion cross section ratios together with the data at several 
energies for ir + (dots) and tt~ (open circles) [||, ^6|. In our NLO calculations 

we used the previously fixed (fcf,) values at scales Q = pr,jet/2 (dashed lines) and 
Q = prjet (solid lines). Here the Cronin parameter is C = 0.35 GeV 2 and v max = 3. 

In the absence of the Cronin effect, these ratios would be identically 1. The 
significant deviation of the data from unity is therefore a clear confirmation of the 
nuclear enhancement in the 2 GeV < pt < 6 GeV transverse momentum window. 
At low pt the ratio is smaller than unity, indicating absorption effects. 

It is clear that the height of the Cronin peak depends on the extra kr broaden- 
ing (quantitatively on the product C ■ h p A(b)), and Fig. ^ tells us that the location 
of the peak depends on the value of (k^) pp . Using our previously fixed parameters, 
NLO calculations at Q = pr.jet (solid lines) give the best fit to the experimen- 
tal data (overlapping with LO curves at Q = pTjet/2, which curves not shown in 
Fig. [| see Ref. [E|). In both cases the values of (k^) pp are very much similar and 
(kx)pp ~ 2 GeV . This characteristic value of {kj) pp is in good agreement with 
the experimental results of Ref. [ p9| . However, NLO calculations at Q —p T j et /2 
(dashed lines) with much smaller (k^) pp (see Fig. ||) give peaks shifted to unphysi- 
cally low values of pr- The sensitivity of the Cronin effect to the scale at the NLO 
level dictates our choice of scale in the present work. 
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3.3. Nucleus-nucleus collisions 

In nucleus-nucleus reactions, where both partons entering the hard collision orig- 
inate in nucleons with additional semi-hard collisions, we do not need additional 
parameters, apply the Cronin enhanced width (Q) for both initial partons. Thus, 

dp J dp 

For CERN energies we follow Fig. || and fix the values of (k^) pp — 1.6 GeV 2 
(E beam = 158 AGeV) and 1.7 GeV 2 (E beam = 200 AGeV) like in Ref. [§ for LO cal- 
culations. We keep Q — prjet and the Cronin parameters v m — 3, C — 0.35 GeV 2 . 

In Fig. the results of our calculation for ir° transverse momentum spectra are 
compared to the WA80 [ |(| and WA98 [ [[l) data on central S + S, S + Au, and Pb+ 
Pb collisions. Top panel shows the spectra, bottom panel displays the data/theory 
ratios. Dotted lines represent the results of the pQCD calculation with nuclear 
effects (shadowing and multiscattering) turned off, while solid lines correspond to 
the full calculation. (To test the effect of shadowing alone, we turned it off; the 
modification resulted in a < 10% downward shift in the spectra.) These calculations 
were performed down to px = 1 GeV; however, since pQCD is a theory of hard 
particle production, it is not expected to describe the data below px ~ 2 GeV. 

In S + S and S + Au collisions the pQCD results reproduce the data at 
Pt > 2.5 — 3 GeV. In Pb + Pb collisions we would expect a similar effect. How- 
ever one can see a definite (upto 40%) overestimate in the calculation. The same 
result was obtained in the LO analysis (see Ref. [ [|). The deviation of the data 
from the pQCD prediction indicates the appearance of a new collective effect in the 
Pb + Pb collision. One candidate for this deviation is "jet quenching" , the induced 
energy loss of high energy jets penetrating hot dense matter. This effect can be 
seen clearly at RHIC energy (see the data in Refs. [ [jffi, |33|, and an explanation 
in Ref. [ |34| ). Since jet quenching is strongly increasing with energy, its moderate 
appearance at CERN SPS in central Pb + Pb collisions is acceptable. 

4. Summary and conclusion 

We presented a pQCD based parton model calculation in next-to-leading order 
powered by a jet-level K factor. The intrinsic transverse momentum distribution 
of the partons inside nucleons was included into the model. The Gaussian width 
of this distribution in pA and AA collisions is controlled by two terms: the free 
pp value {k T )p P , fixed by several different experiments, and a nuclear part, which 
gives extra enhancement due to semihard collisions. (Note that at NLO level the 
factorization scale was chosen to be Q = pT,jet-) While the nuclear part controls 
the height of the Cronin peak, the {k^) pp part determines its position. 

We analyzed CERN SPS data and found that in Pb + Pb collisions the ex- 
perimental data are below the pQCD prediction by 30 — 40 %. This is similar to 
what was seen in LO calculations [ |). With the present choice of scale and NLO 
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Fig. 7. Neutral pion production compared to data from SPS experiments WA80 
and WA98 for central heavy ion collisions. Top panel shows the data and the calcu- 
lated NLO invariant cross sections with the nuclear effects turned off (dotted lines) 
and for the full calculation (solid lines). Bottom panel displays the data/theory 
ratios for the full calculation. 

approximation, the overestimation appears to be independent of the order of the 
pQCD calculation and indicates some collective nuclear effect. One candidate is 
the induced jet energy loss in hot dense matter ("jet-quenching"), which becomes 
significant in central Pb + Pb collisions. 
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